Stability of Excited Dressed States with Spin-Orbit Coupling 
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We study the decay behaviors of ultracold atoms in metastable states with spin-orbit (SO) cou- 
pling, and demonstrate that there are two SO-coupling-induced decay mechanisms. One arises from 
the trapping potential and the other is due to interatomic collision. We present general schemes 
for calculating decay rates from these two mechanisms, and illustrate how the decay rates can be 
controlled by experimental parameters. We experimentally measure the decay rates over a broad 
parameter region, and the results agree well with theoretical calculations. This work provides an 
insight for both quantum simulation involving metastable dressed states and studies on few-body 
problems with SO coupling. 
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Recently, synthetic magnetic field and a restricted class 
of spin-orbit (SO) coupling have been successfully real- 
ized in ultracold Bose gases [lrfl| and degenerate Fermi 
gases 0, Q . In addition, many schemes have been pro- 
posed to create general gauge fields, particularly the 
isotropic Rashba SO coupling; see Ref. M for a review. 
These will bring about novel quantum systems of SO 
coupled atoms that display many interesting phases fiol - 
16[, and provide a new and powerful tool for quantum 
simulation. However, in many of these proposals, one 
or more metastable states, e.g. dark states, play essen- 
tial roles. Thus, the lifetime of atoms in the metastable 
states becomes crucial in practical implementation of 
these schemes. 

In this Letter, we carry out a thorough and quantita- 
tive study of the decay behavior of excited dressed states 
with SO coupling. Wc find that, due to the SO coupling, 
decay mechanisms can arise both from single-atom mo- 
tion in inhomogeneous trap potential and from two-body 
collisions. The trap-induced decay rate is determined by 
the trap frequency, while the collisional decay rate is con- 
trolled by atomic density and scattering length. Further, 
we present rigorous methods for calculating decay rate 
for each mechanism. Our method for the rate of colli- 
sional decay can be applied to both elastic and inelastic 
scattering of two atoms with a general type of SO cou- 
pling and arbitrary scattering length. Finally, we experi- 
mentally investigate the decay behavior of a SO coupled 
87 Rb BEC prepared in a metastable state for a broad pa- 
rameter region. The experimental results are compared 
to theoretical calculations, and the excellent agreement 
supports the validity of our theory. 

Our work provides a comprehensive understanding of 
the SO-coupling-induced decay and thus can serve as a 
valuable reference for experimental realization of propos- 
als involving metastable dressed states. For a given sys- 
tem of interest, one can use our theory to figure out which 



FIG. 1: (a and b): Sketch of the trap-induced decay mecha- 
nism (a) and the collisional decay mechanism (b). (a) illus- 
trates two steps of the trap-induced decay, and (b) illustrates 
the two typical processes in the collisional decay. 



is the dominant mechanism, and then apply appropriate 
approaches to control the stability of those metastable 
states. In addition, our method for the exact calcula- 
tion of interatomic scattering amplitude sheds lights on 
few-body problems with SO coupling. 

The trap-induced decay. The Hamiltonian of a single 
atom with SO coupling takes form Hi^ = p 2 /(2m) + 
M(p) + V(r) = Hq(p) + V(r), with m the atomic mass, 
p and f the atomic momentum and position operator, re- 
spectively, and V(r ) the trap potential. The SO coupling 
is described by the p-dependent operator M(p). For in- 
stance, for effective spin- 1 /2 systems in Refs. [3-[! 0, S , 
one has M(p) = Sa z /2 + Q.a x /2 + 2k r p x a z , where a is the 
Pauli operators, S is the two-photon detuning, k r is the 
recoil momentum and f2 is the Raman-coupling strength. 

It is clear that the eigen-state of Hq is |k)|a(k)), where 
|k) satisfies p|k) = k|k) and the state |ct(k)) in the spin 
space is the eigen-state of M(k). If there were no SO 
coupling, both M(k) and |a(k)) are k- independent. In 
this case, since the trap potential V is independent of 
the atomic spin, it cannot induce the transition between 
two eigen-states with different a, or the decay from the 
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excited spin state. In the presence of SO coupling, both 
M(k) and |a(k)) depend on k. Thus, |a(k)) and |a'(k')) 
with k ^ k' can overlap with each other even if a a'. 
Due to this fact, V(r ) will couple two dressed states with 
either different a, or different k, or both. Therefore, 
excited dressed states can decay to lower-energy branch 
via the single-atom motion in trap V. 

Here the trap- induced decay process can be understood 
as two steps as illustrated in Fig. 1(a). First atoms tunnel 
from the initial state to the energy-conserved states in the 
lower branch (the solid arrow in Fig. 1(a)). The rate Tib 
of this process can be calculated by Fermi's golden rule 
(FGR). Second, due to the dissipation effects given by 
the environment, the atoms further decay to states with 
lower energy (the dashed arrow in Fig. 1(a)). Usually, 
the second process is much faster than the first, and thus 
the total rate is given by Tib- 

Now we calculate Tib in the momentum representa- 
tion. The atomic state \ip{t)) at time t is described by 
the spinor wave function \ip(k,t)) = (k.\ip(t)), and the 
position operator f is proportional to the gradient of the 
momentum k, i.e., f = iVk- The harmonic trap poten- 
tial can be written as V = —J2j= x y z ( m /2)uj 2 d 2 /dk 2 , 
which behaves as the "kinetic energy" of the atom motion 
in the k-space. Furthermore, \ip(k,t)) can be expressed 
as \ip(k,t)) = Y^ a ^(k, £)|a(k)). Then the Schrodingcr 
equation id\ip(k,t))/dt = [V + Ho(k)]\tp(k,t)) can be re- 
written as the equations for each component -0 Q (k, t): 
idipa/dt = J2pT a piJjp + £ a Qs.)il> a , where £ Q (k) is the 

eigen-energy of Hq for |k)|a(k)), and T Q( g is defined 
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i6 a pd/dkj+i(a(k)\d/dkj |/3(k)). It is clear that the terms 
i(a(k)\d / dkj\f3(k)) play the similar role as the effective 
gauge field in the Born-Oppenheimer approximation fl7j |. 
The decay of the atoms from the excited dressed state, 
or the transition between dressed states with different a, 
is induced by these terms with a^/). 

Suppose that the initial wave function of the atom 
is |V>i(k)) = 0(k)|a(k)), with 0(k) satisfying [T aa + 
£„(k)]0(k) = £<£(k). Then r lb can be given by FGR 
as 



r lb = 2ttJ2pp( £ ) 



dk<i>*Jk)Tp a cj>(k) 



(i) 



where ^(k) satisfies \Tpp + f l g(k)](/> | g(k) = ecppik) and 
pp(e) is the associated density of states. Since V oc 
ui^. yz , the decay rate Tib is proportional to z , and 
thus can be controlled by tuning the trap frequency. 

The collisional decay. For two SO-coupled atoms un- 
der scattering, the general Hamiltonian takes form H 2 \> = 
ffo(l) + H (2) + U(f 12 ) = H F + U, with ff (i) for the 
free motion of the ith atom (i = 1,2) and U(yi 2 ) the 
interaction potential of the two atoms. Here r\ 2 is the 
relative position operator of the two atoms. We shall 
consider the simple case where U is independent of the 



atomic spin. If there were no SO coupling, U cannot 
induce transition between different spin states or the 
atomic decay from excited spin states. In the presence of 
the SO coupling, the free-motion state of the two atoms, 
or the eigen-state of Hf , becomes the two- atom dressed 
state |c) = |ki)i|ai(ki))i|k 2 )2|a2(k2))2- Here we define 
c = (ki, ai, k2, a 2 ) as the set of the four quantum num- 
bers. Similarly as in our above discussion, states |a(k)) 
and |a'(k')) can overlap with each other when a =/= a' . 
Then we have (c\U\c') ^ even if (01,02) ^ (ai,a' 2 ), 
and U can introduce inelastic collisions or the transi- 
tions between the states with different quantum number 
(011,0.2). This leads to the decay of atoms from excited 
dressed states (Fig. 1( b)). 

The above discussions are applicable to both bosonic 
and fcrmionic systems. Hereby we shall focus on a system 
of bosonic atoms condensed in an initial dressed state 
|k )|ao(k )) with atomic density n . The characteristic 
decay rate T2 b for the collisional decay is defined as T2 b = 
noK. Parameter K is the inelastic-collision rate defined 
as K = 2crv, where a is the total cross-section of the 
inelastic collision, v is the relative velocity of the two 
atoms before collision, and the factor 2 comes from the 
bosonic statistics. According to the standard scattering 



theory [18(, the factor K is given by 
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(o4, 04)^(0*1, a 2 ) 



dkidk^f 5 K \f(c',co)\ 2 (2) 



with Co = (ko, cxq, ko, ao) and /(c',c) is the scattering 
amplitude between the incident state \c) and the output 
state |c'} with d = (k[, a[, k' 2 , a' 2 ). The Dirac S function 
5g means that the total energy is conserved during the 
scattering process, and <5k = S(k± + k 2 — k[ — k' 2 ) is for 
the momentum conservation. 

In Eq. d2]), the scattering amplitude / is defined as 
/(c',c)<5k = —2ir 2 m{d\U\c+), with |c+) the scattering 
state given by the Lippmman-Schwinger equation (LSE) 
\c+) = |c) + GoU\c+). Here the free Green's operator Go 
is defined as G = [£ ai (ki) +£ a2 (k 2 ) + i0+ - Hp]' 1 . We 
assume U(ri 2 ) is a short-range potential and becomes 
negligible when |ri2 1 is larger than an effective range 
In this case, when there were no SO coupling, the wave 
function (i"i2|c+) of the scattering state should satisfy 
the Bethe-Periels boundary condition (BPC) (ri 2 |c+) oc 
(l/|ri2 — 1/a) in the short-range region i?» < | r- 12 1 << 
1/fc, where k = |ki — k 2 |/2, a is the scattering length 
and |ri 2 ) given by fi2|ri 2 ) = ri 2 |ri 2 ). Therefore, one 
can first derive the expression of (ri2|c+) in the short- 
range region with the LSE, and then substitute the result 
into the BPC and obtain l^) as well as /. 

In the presence of SO coupling, the Hamiltonian of two 
colliding atoms is revised in the whole range of the inter- 
atomic distance |ri2 1 , and thus the few-body properties 



are also strongly affected by SO coupling [6|, 19|-23|. In 
particular, the Bethe-Periels boundary condition needs 
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FIG. 2: Decay behavior of the BEC of 87 Rb atoms m our ex- 
periment. The measured fraction lZ(t) of condensate versus 
different hold time for cases with Q = 0.6E T , S — 6£" r (a) 
and SI = 0.9i? r , S = — 6E T (b), with the dispersion relation- 
ships shown in the insets of (a) and (b), respectively. The 
blue curve is obtained by fitting the experimental data with 
our theoretical function of lZ(t). 



to be modified [24J and the free Green's function will be 
re-calculated. Nevertheless, the above procedure is still 
valid. With this, we find that f(c', c) is given by [24 1 



/= 2(4(^)1^(^)1- 
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FIG. 3: The characteristic time r = l/(Kn c ) of the collisional 
decay in our experiments. (a,c): r as a function of fl for 
S = 6-E r (a) and 8 — —6E r (c). (b,d): r as a function of S 
for 5 > (b) and S < (d), with fi = 1.2E T . The values of 
r obtained from experiments (open circle with error bar) are 
compared with the theoretical calculation with Eqs. Q and 
© (blue solid line). 
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of SO coupling, 
pressed as T = im/ (Aix)^/E — (ki + k 2 ) 2 / (4m) — 
(2vr)- 3 E A)fe I dp (A" 1 - A" 1 ) \[ J=l 2 \^)W 3 {^)\, 
with E = £ ai (ki) + f Q2 (k 2 ), Pli2 = (ki + k 2 )/2 ± p, 
A = E + i0+ - %(pi) - £/? 2 (p 2 ) and A = 
£ + i0+ - (pi + p!)/(2m). When m/(47ra) 
is much larger than eigenvalues of J-, one has 
/ = -a(a 2 (k 2 )|a 2 (k 2 ))(a' 1 (ki)|ai(ki))-i.e., the re- 
sult given by the FGR approximation For large 
a, contribution from J- becomes significant, and the 
FGR fails. With the expression ((3|) of the scattering 
amplitude, one can obtain the factor K and the decay 
rate T 2 b = UqK. 

Experiment. Our experimental layout has been de- 
scribed in Ref. 0. A BEC of 2.5 x 10 5 87 Rb atoms in 
the F = 1 manifold is created in an optical dipole trap 
with frequencies of {oo Xl us yi lj z } = 2n x {30, 30, 50}Hz. 
The single-atom Hamiltonian in the x-direction is given 
by H x = H 0x + V, with 
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and V = mujlx 2 /2 the dipole trap potential. Here 11 
is the strength of Raman coupling, 5 is the two-photon 
Raman detuning and e is the quadratic Zeeman shift 
given by a homogeneous bias magnetic field. Symbol k r 
represents the recoil momentum, and E r = fc 2 / (2m) = 



2ir x 2.21kHz is the recoil energy. Diagonalization of the 
Hamiltonian Hq x leads to three momentum-dependent 
eigen-states \k x )\a(k x )) (a = 0,±1) with eigen-energies 
£-i (k x ) < So (k x ) < £+i (k x ). Two examples of the dis- 
persion curves arc shown in the insets of Fig. 2(a,b). 

For experiments with 8 > 0, the BEC is first prepared 
in the bare state \F = l,m_p = —1), and transformed to 
the \m,F = 0) state with a 7r-pulse. Then we adiabatically 
turn on the SO coupling, so that the BEC is prepared 
in the middle dressed state with a = and k x around 
some value ko, corresponding to the global minima of the 
£o (k x ) curve for S > (the inset of Fig. 2(a)). In the 
experiments with S < 0, the SO coupling is adiabatically 
applied on the BEC in the state \F = l,mp = —1), and 
then the system is prepared in the dressed state with 
a = and k x around the local minima for S < (the 
inset of Fig. 2(b)). The Raman coupling is held for a 
variable duration th- During this time interval, the atoms 
can decay from the initial dressed state with a = to 
those with a = — 1 0, Q. At t = th_, the Raman lasers 
and the dipole trap are suddenly turned off. With the 
Stern-Gerlach technique, a time-of-flight image is taken 
to measure the number of atoms remained in the BEC 
after the decay process. As an example, the fraction of 
remaining atoms is shown in Fig. [2] as a function of th- 

Data Analysis. The numerical calculations with Eqs. 
(|1I3|) show that in our system the rate Tib and T 2 b of 
the trap-induced and collisional decay are of the order 
10 _5 Hz and 10Hz, respectively. Therefore, the trap- 
induced decay in our experiments is negligible, and we 



4 



150r 



.§ 100 



B 50 



(a) 



■ 1/r „ 
■1/T„ 



40 



50 60 70 80 
(o/2n (Hz) 




FIG. 4: The characteristic times l/Tib (blue solid line) and 
l/r2b (magenta dashed line) of the trap-induced and colli- 
sional decay of the system in current experiments. The pa- 
rameters are set as £7 = 0.2E T and S = 4E r , and the conden- 
sate density is taken as 2.9 x 10 13 cm _,i . We plot the decay 
times as a function of the trap frequency u x with the scat- 
tering length a = 100ao (a), and a function of the scattering 
length a with u) x = (27r)60Hz (b). 



can safely consider the inelastic scattering only. When 
two ultracold atoms with a = decay to the a = — 1 
branch via scattering, they likely become thermal due 
to the large energy gap, which however, is not suffi- 
cient for atoms to escape from the trap. These ther- 
mal atoms will also collide with condensed ones. Thus, 
the decreasing of the ultracold-atom density n c can be 
described by dn c /dt = —Kn\ — Ln c (no — n c ), where 
no ~ 2.9 x 10 13 cm~ 3 is the initiate atomic density of the 
BEC, K is defined as before and the parameter L repre- 
sents the collision rate between a thermal and an atom 
in the condensate. Then the BEC fraction lZ(t) in the 
atomic cloud can be obtained from that equation and ex- 
pressed as H{t) = L/[L + K{e tn ° L - 1)] . According to 
this expression, the characteristic decay time of lZ(t) is 
r = 1/T2b = l/(n>oK) when L < K. The experimental 
data, as shown in Fig.[5J were fitted by this function with 
both K and L as fitting parameters. Our fitting result 
shows that the condition L < K is satisfied in our exper- 
iments. The lifetime r obtained from our experimental 
results are shown in Fig. [4] We further theoretically cal- 
culate the coefficient K with Eqs. ([2]) and ([3]). The calcu- 
lated values of the lifetime r are also plotted in Fig. 0] as 
blue curves. The agreement between the theoretical and 
experimental results is excellent |25j . This confirms our 
analysis of the decay mechanisms and the calculations of 
scattering amplitude. 

Discussion on the control of stability. In this letter we 
have shown the two mechanisms of decay in the ultra- 
cold gases with SO coupling, carried out the calculation 
of two decay rates, and presented a comparison with ex- 
periments. This guides us how to control the stability 
of excited dressed state in current setup. For instance, 
as shown in Fig. 4, we plot characteristic times 1/ Tib 
and I / T2b of the trap-induced and collisional decay as 
functions uj x and a for the SO coupling realized in cur- 
rent experiment. When ui x > (2-7r)60Hz, as shown in Fig 



4 (a) , the decay is dominated by the trap- induced decay, 
and then the decay rate can be controlled by the trap fre- 
quency As the trap frequency decreases, the lifetime of 
excited state gets longer. However, when lu x < (2tt)60Hz, 
the collisional decay becomes dominating. In this region, 
the decay rate is no longer sensitive to the trap frequency, 
but it can be controlled by the atomic density and the 
scattering length. For instance, as shown in Fig 4(b), 
when lu x = (27r)50Hz, the life time determined by the 
collisional decay can be increased by reducing the scat- 
tering length a or atomic density. 

As emphasized before, our analysis and calculation 
method are very general and can be applied to ultra- 
cold gases with any kind of SO coupling. Thus, similar 
stability analysis as discussed above for current exper- 
imental system can be straightforwardly carried out for 
other realizations of SO couplings. It will help the imple- 
mentation of those schemes involving metastable states, 
and therefore open up more possibilities for rich physics 
with synthetic gauge field. 
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